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ABSTRACT 

66" 

Q\ , In this Letter we explore a version of the test of cosmological geometry pro- 

posed by Alcock & Paczyhski ( |1979| ), using observations of the Lyman-a forest 
in the spectra of close quasar pairs. By comparing the correlations in absorp- 
tion in one quasar spectrum with correlations between the spectra of neighboring 
quasars one can determine the relation of the redshift distance scale to the angle 
distance scale at the redshift of the absorbers, z ~ 2 — 4. Since this relationship 
depends on the parameters of the cosmological model, these parameters may be 
Q\ [ determined using the Lyman-a forest. While this test is relatively insensitive to 

1^ the density parameter VL m in a dust-dominated universe, it is more sensitive to the 

O \ presence of a matter component with large negative pressure (such as a cosmo- 

logical constant A) and its equation of state. With only 25 pairs of quasar spectra 
"^h at angular separations 0.5' — 2', one can discriminate between an Vt m = 0.3 open 

universe (A = 0) and an Q m = 0.3 flat (A-dominated) universe at the 4 — a level. 
The S/N can be enhanced by considering quasar pairs at smaller angular separa- 
tions, but requires proper modeling of nonlinear redshift space distortions. Here 

ci' 



the correlations and redshift space distortions are modeled using linear theory. 



Subject headings: cosmology: theory — intergalactic medium — quasars: ab- 
sorption lines — large-scale structure of universe 



1. Introduction 

Recent supernova la observations have generated a lot of interest in cosmological mod- 
els where a significant fraction of the energy contents has negative pressure ( [Ferlmutter 
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et al. 1997; Kiess et al. 1998 ; l^arnavich et al. 1998|) . A common way of parameterizing 



the equation of state of this component, which we call Q, is p = wp, where p and p denote 
the pressure and density respectively. The cosmological constant A corresponds to the case 
w = —1. 

It is important to have independent ways to constrain the abundance and properties of 
such a component, as different methods suffer from different systematic errors, and, perhaps 
more importantly, different methods are sensitive to different combinations of parameters. 
In this paper, we discuss a version of a test proposed by Alcock & Paczyhski (|1979fc AP 
hereafter), which is particularly sensitive to the presence/absence of Q. They observed that 
an object placed at a cosmological distance would have a definite relationship between its 
angular and redshift extents, which is cosmology-dependent. 

Consider an object with mean redshift z, and angular size 9 which thus has transverse 
extent (in velocity units) 

u ^ = rh DA ^ e ■ w 

Here H is the Hubble parameter at redshift z, and Da(z) is the angular diameter distance. 
For spherical objects the radial and transverse extents are equal, but more generally if the 
object is squashed radially by a factor a s , then the radial extent is u\\ = ^§ = a s u± Here 
c is the speed of light and u±, u» <C c is assumed. A plot of u±/9 is shown in Fig. |[ Note 
how its value for a Q-dominated universe (for w H — 1/3) differs significantly from that for 
a no-Q-universe. 

Various incarnations of this test have been discussed in the context of galaxy and quasar 
surveys (e.g. [Ryden 19951 ; |Ballinger et al. 19961 ; PVIatsubara fc Suto 1996| ; |Popowski et al. 1998 



De Laix fc Starkman 1998| ) where the "object" used is the two-point correlation function, 



whose "shape" (a s above) need not be spherical because of redshift-anisotropy induced by 
peculiar motion. In the case of the Lyman-alpha (Lya) forest, we cannot observe the full 
three-dimensional (3D) correlation directly. Instead we can measure the one- dimensional 
(ID) correlation along a line of sight (LOS), and the cross-correlation between two (or 
multiple) close-by LOS, or their Fourier counterparts: the auto- and the cross-spectra. The 
two are related to each other through an underlying 3D power spectrum. These relations 
are spelled out in §§ (for auto-spectrum, see also [Kaiser fc Peacock 1991| ). Using a method 



developed by Hui ( |1998a|) building on earlier work by Croft et al. (|1998|) , one can invert the 



auto-spectrum to obtain the underlying 3D power spectrum, which then allows one to predict 
what the cross-spectrum should be, except that the prediction is cosmology-dependent. A 
comparison between the observed and the predicted cross-spectrum is then our version of 
the AP test. It is not the only version possible in the context of the forest, but we will use 
this version to gain some intuition on the sensitivity of the AP test to various cosmological 
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parameters. A hypothetical implementation of this test is studied in §3.2. 

An important problem in this application of the AP test is the modeling of the redshift- 



space distortions (the "shape"). We give a first estimate using linear theory in § [3.1j We 
conclude in §|4] with an assessment of the expected S/N of the AP-effect-measurement in the 
Lya forest, and remarks on lines of further investigation. 



2. Auto- and Cross-spectrum: Definitions and General Formula 

Given the observed transmission / = e~ T as a function of the velocity u\\ along a 
LOS and the angular position 9 on the sky, let us define 5f = (/ — /)// where / is the 
mean transmission. Then the auto-correlation along the LOS and the cross-correlation 
!;(. between two close-by LOS of angular separation 9, and their Fourier counterparts the 
auto-spectrum P,f and the cross-spectrum P^, are respectively defined by 

(5 f (u\ l ,6>)5 f (u\ l +u lb 6'))=t{ l (u ll ) , (5 f (u\ l ,6')5 f (u\ l +u h 6' + 6))=d(u lb 6) (2) 

P((k { \) = J £f(«||)e- ifc 'i"M M || , Pf(k h 9) = J d(u h e)e- ih nd Ull (3) 

where we use fen to denote the wave- vector along the line of sight. 

The two quantities are related to an underlying 3D power spectrum pf(kn,k), where 
we use ~ to denote the 3D nature, and k is the length of the 3D wave- vector: 

p{(h) = I P f (h^)k^ (4) 
pl(k h e) = / p f (k h k)j [k ± u ± {e)]k— 

where Jq[x] = Jq W e~ txcosa da/ (2tt) is the zeroth order Bessel function, k± = ^Jk 2 — kp and 
u± is the transverse velocity-separation for the given 9 (eq. [HI). 

We deliberately allow P*{k\\,k) two arguments to account for the possibility of an 
anisotropic 3D power spectrum, as for instance in the presence of distortions by peculiar 
motion, or thermal broadening (bear in mind that line-broadening always acts along, not 
transverse, to the LOS). The tricky part is how to model this distortion. A general form is 

P f (k h k) = W(k\\/k,k)P f (k) (5) 

where W is a distortion kernel, and P-'(k) with only one argument represents the transmission 
power spectrum in the absence of peculiar motion and line-broadening. 
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As we will discuss in §3^, eq. (pj) together with eq. (ffl) allow us to predict the cross- 
spectrum, once the auto-spectrum is given and W as well as the relevant cosmological pa- 
rameters determining the angle- velocity relation (eq. |]T|) are specified. A comparison of 
the predicted and observed cross-spectra then allows one to discriminate between different 
cosmological models, which is our version of the AP test. We will perform a first estimate 
of W using perturbation theory. 



3. A Perturbative Estimate 

3.1. Perturbative Auto- and Cross-spectra 

A linear perturbative calculation (i.e. assume 5f <C 1) of the auto-spectrum has been 
carried out in Hui ( |1998a|) . The cross-spectrum can be computed in a very similar manner. 
We simply state the results here, and refer the reader to Hui (|1998a ) for details. 



Essentially, both power spectra follow from the substitution of eq. ([|) into eq. (|3]), with 
W given byfl 

\ 



W(kJk,k) 



1 + ^ + 



exp[-k'i/kf] (6) 



with 



Pf ~ 2 - 0.7( 7 - 1) ' Ab ~ 8-2.8(7-1)^ (?) 

where fa = dlnD/dlna with D being the linear growth factor and a the expansion scale 
factor, bT is the thermal broadening width associated with a mean temperature of T , and 
7 specifies the temperature-density relation through T oc p 7_1 where T is gas pressure and 
p is the density. The smoothing scale k» is associated with the thermal broadening scale as 



well as observation resolution, (see Hui |1998a| for details). 



An important feature of eq. (||) is that on large scales (k\\, k k», 1/&t ), it reduces to 
the famous Kaiser ( |1987| ) formula, with (3f playing the role of the galaxy-bias factor: 



2\2 



W(k\\/k, k) -> W L s(k {l /k, k) = (1 + Pfk^/k 2 ) 



As shown in Hui (|1998a ), this allows a one-parameter-only (/?/) inversion of the large 



scale 3D isotropic power spectrum pf(k) using this linear integral equation, which follows 



3 We have equated (k) with the A'P p (k) exp[— k 2 /kp] in Hui ( 1998a| ), where A' is a constant, fc_F is 



roughly the Jean scale, and P p (k) is the 3D real-space mass power spectrum. 
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from eq. and (|): 

lf(*||)= / k)P f (k) — — . (9) 

The above, in discretized form, can be viewed as a matrix equation, and one simply inverts 
a matrix proportional to W (which is in essence upper or lower triangular because of the 
limits of integration/summation) to obtain P*(k) from the auto-spectrum. It can be shown 
that Wis (eq. ||) could be used instead of the full W in obtaining the large scale P^(k), 



with good accuracy. The reader is referred to Hui ( 1998a ) for illustrations of this method. 

An interesting bonus of the distortion kernel given in eq. @ is that the factor of 
exp[— /c|//cfj 2 ] is commonly used to describe nonlinear redshift distortions of galaxy distribu- 
tions (see e.g. Peacock fc Dodds 1994|) . This should come as no surprise because nonlinear 
distortions arise from virialized objects which have a Maxwellian distribution of velocities - 
precisely the form for thermal broadening as well, which determines fcjj. One can then use 
the perturbative kernel in eq. (|6]) but allow fcjj to vary to account for nonlinear distortions. 
We will not pursue that here, but will come back to it in §[|. It should be emphasized that 
even if W turns out to be substantially different from eq. (j^) or (|j), perhaps due to nonlin- 
ear clustering, the procedure of inverting a triangular matrix using eq. (||) should generally 
work. 

Given the 3D isotropic power spectrum pf(k), one can predict the cross-spectrum using: 

/*oo _ hrlk 

P{{k h 6)= W (k {l /k,k)pf(k)J [k ± u ± (6)}— (10) 

which follows from eq. (|4]) and (|5|). As in the inversion procedure above, we will use Wis 
instead of the full W, when performing the above integration. We will see in the next section 
that the induced error in the predicted large scale cross-spectrum is small. Note that the 
rapid oscillation of Jo an d decay of P^ at high kj_ or k means the large scale P£ is not 
sensitive to assumptions made about the small scale distortion kernel. 



3.2. The Alcock-Paczyriski Test 

The AP test for the Lya forest could be implemented in several different ways. The 
version adopted here should be seen as a first step towards understanding the sensitivity of 
the test. We will discuss the issue of other possible implementations in §|4]. 

Suppose one is given a set of idealized observed transmission power spectra, say the 
auto-spectrum and the cross-spectra for two different angular separations 6 = 1' and 2'. In 
real life, one of course does not know a priori the underlying cosmological model behind this 
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set of observations. For the sake of our testing here, let us construct these power spectra 
by assuming an input "true" model: the open Cold-Dark-Matter (OCDM) universe with 
fl m = 0.3 and no Q. The full distortion kernel (eq. 0) is used to compute these power 
spectra, with 7 = 1.5, and &t , kf< and hp corresponding to a gas of temperature 10 4 K ( see 



Hui |1998a| for details). The mass power spectrum is assumed to have a shape parameter of 



r = 0.25. Let z = 3. 

Pretending we have no insider information on the underlying cosmology, our version 
of the AP test comes in two steps. First, perform the inversion from the "observed" auto- 
spectrum P,| (k\\) to the 3D isotropic power spectrum pf(k) using eq. (|9|) and Wls i n ec i- 

By using the latter instead of the full kernel as in eq. (^|), we will be making an error. 
However, the error will be small on sufficiently large scales, as we will see. The second step 
involves predicting the cross-spectra for the corresponding angular separations, using eq. 
(10), once the 3D isotropic power spectrum P*{h) is computed. 



Both steps require the assumption of a cosmological model. In the first step, the cosmo- 
logical density parameters determine the amount of redshift-space anisotropy that needs to 
be accounted for, through the parameter (3 f or fa (eq. || or 0). The temperature-density- 
relation index 7 also enters into /3/; we will come back to this below. In the second step, 
cosmology dictates the angle- velocity relation through the quantity u±(9) (eq. Jl| & [|Kj). 



To parameterize the cosmological model that we need to assume in carrying out the 
AP test, we use Q m , Qk and Qq, which correspond to the mass, curvature and Q density 
parameters. They sum to 1. The Q component is described by an equation of state of 
the form p = wp where — 1 < w < 0. The limiting cases of —1 and correspond to the 
cosmological constant A and matter m respectively. For w > — 1, the Q component could 
cluster (e.g. |Frieman et al. 199"5] ; [b'erreira fc Joyce 1998| ; |Caldwell et al. 19"9~T| ), which would 



affect the distortion parameter We consider the simple case of no Q-clustering here (see 
e.g. Turner fc White 19971 ). 



Assuming any set of parameters different from that of the underlying input model 
(OCDM) will result in predicted cross-spectra different from the "observed" ones. The 
fractional errors, for a set of 5 cosmological models, are shown in Fig. [I]. 

In the case where the actual input cosmology is used, the cross-spectra are quite accu- 
rately predicted. It is not perfect, especially at smaller scales, because the inversion proce- 
dure in the first step of our version of the AP test is inherently approximate by avoiding the 
modeling of small scale distortions (see |Hui 1998a|) . However, in cases where the "wrong" 
model-cosmology is assumed, the predicted cross-spectra are systematically different from 
the "observed" or input cross-spectra. 
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The effect is most pronounced when Q is present. The canonical example is the cosmo- 
logical constant with w = —1, which gives the strongest departure, among the five models, 
of the predicted cross-spectra from those of the input cosmology, with differences as large as 
200% at sufficiently large k\\'s. Increasing w while keeping Qq fixed (the w = —1/3 model), 
or decreasing Qq while keeping w fixed (adding in mixture of curvature, as in the somewhat 
perverse model of Q m = 0.2, Q k = 0.4 Sz Q\ = 0.4) tends to bring the predictions closer to 
that of the input open universe, but still with substantial differences. The Q m = 1 critical 
matter density universe is closest to the input model (that is, aside from the input model 
itself), with a difference of about 10 — 20%. 

As we discuss before, the cosmological model influences the prediction of the cross- 
spectra through two parameters: /n (or /3f] see eq. |7|]) of the distortion kernel, and u±(9) 
(eq. [!]]) of the velocity- distance relation. To understand approximately the size of the effect 
we are measuring, it is simplest to ignore redshift-distortion first, and consider the influence 
of the second parameter alone. 

Putting W = 1 and assuming a power- law P*(k) = Bk n in eq. (|T0"D, it can be shown 
that ( [Abramowitz fc Stegun 1964|) 



H(k\\,9) = ^ [2*||/«l] ' Ki + ±(h\\u L )IT(-nl2) (11) 
where K u is the z/-th order modified Bessel function and V is the gamma function. K u (x) has 



an asymptote of y7r/(2x) exp[— x] in the large-x limit. ( For v = ±1/2 this is exact.) For 
the scales of interest, n should vary somewhere between —2 and —3. For sufficiently large 
fey (i.e. k\\u± ^7r), the u± and h\\ dependence of the cross-spectrum is given by: for n = —2, 
Px(k\\,9) oc exp[— k\\u±]/ Jk\\u±, while for n = —3, P£(k\\,9) oc exp[ — Aj||itj_]/A;||. 



This exponential dependence on u±(9) explains the large differences between different 
models at the smaller scales. (Of course, the exponential suppression of power at high fcy's 
also means the cross-spectrum is going to be hard to measure at such scales. See §f|.) We 
show in Fig. ^| two curves depicting how u±/9 varies with w and with Q^. It can be seen 
that pure matter-plus-curvature models (with no Q) generally have a higher u±/9 than 
Q-dominated models. Combining such information with Fig. |TJ, we can say that the cross- 
spectra of models with a significant Q component are generally 50% higher, or more, than 
those of no-Q-models, at k\\u± ~ tt. This can be understood using our analytic formula in eq. 
(p|), assuming n ~ — 2 to —3. 

Put in another way, given an observed auto-spectrum, the inferred correlation length 
in physical distance scale is larger in a universe dominated by Q compared to a no-Q- 
universe; ignoring redshift-space distortions, this implies a larger cross-correlation length, 
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and manifests itself in a stronger cross-correlation for the Q-dominated model. That is what 
we see in Fig. |T[ 

A good illustration of redshift-space anisotropy can be seen by a comparison of the 
fl m = 1 model and the open model with Q m = 0.3 & Qk = 0.7. According to the dotted line 
of Fig. [|, the two have very similar u±s with the latter's a little smaller. This means, with 
no redshift distortions, the latter should have a cross-spectrum close to, but above, that of 
the Q m = 1 model. Exactly the opposite is observed in Fig. [3]. This is because peculiar 
motion induces a stronger anisotropy of the correlation function in the Q m = 1 model, in 
the sense of a greater squashing of the correlation length along the LOS (see eq. ||). The 
end-result is a stronger cross-correlation of the critical-matter-density model over the open 
model. The effect is not large. It shifts the cross-spectrum of the Q m = 1 model relative to 
that of the Q m = 0.3 - Qk = 0.7 model by about 10%. A similar conclusion holds for most 
other models in that redshift-space distortions change their fractional differences by about 
10 — 20%. For reference, the values of fn (eq. [0]) for the 5 models in Fig. [3] are 0.98, 0.77, 
0.78, 1.0 and 0.77 from top to bottom . 

It can also be seen from our hypothetical AP test that besides constraining the absence 
or presence of Qq, the predicted cross-spectra can also be used to discriminate between 
different equations of state in cases where Qq is known. Unfortunately, according to Fig. ^ 
the quantity uj_/6 takes rather similar values for w < — 1/3, which are also values that seem 
to be allowed by current supernova observations ( |Garnavich et al. 1998| ). The test does have 
good discriminating power among models with larger u>'s, however. 

The distortion parameter /?/ depends on, in addition to the various JTs, the temperature- 
density-relation index 7 (eq. JF|), which while not known precisely, has been shown by Hui & 
Gnedin ( |1997| ) to generally lie in the range 1.3 < 7 < 1.6. Varying 7 within this range changes 
the cross-spectra by < 10%. 



4. Discussion 

At this point, a natural question to ask is: how well can we measure the cross-spectrum? 
And given the expected signal-to-noise of such a measurement, what is the expected discrim- 
inating power among different cosmologies? 

The expected variance in the measured cross-spectrum for a single /cy-mode is approx- 
imately given by (P<f ) 2 + (P*) 2 (see [Hui 1998~5| ). This is obtained by ignoring shot-noise, 



which is likely to be accurate because the Poisson- variance is known to be sub-dominant com- 
pared to sample- variance (at least for the brighter quasars, see e.g. |Uroft et al. 1998| ), and 
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assuming Gaussianity, which is probably a good approximation by the central-limit-theorem, 
especially on large scales ( [Hui 1998b[ ). Since Pi 3> P^ in general, we can approximate the 
signal-to-noise for each wave-mode of the measured cross-spectrum by P^/Pu- 

On the other hand, our calculation in the last section tells us what 5P^/P{ (Fig. |I|) 
is between any pair of cosmological models. Let us denote these two models by A and B. 
As argued in § |3.2| , the dominant contribution to the difference in Px' s can be estimated by 
ignoring redshift-space distortions. In that case, the difference in the cross-spectra is due 
entirely to the difference in the transverse velocity-separations (eq. U), let us call them uj_ 
and u*l respectively. Assuming n = —3, eq. (|TT1) & (P) allows us to compute both 5Pf</P{ 



and P£/Pl, from which we can deduce the following: one can rule out model B, if A is the 
true model, with a cx-level or signal-to- noise {S/N) of 

S 2 

- = V [exp[-A;||K - <)] - l] 2 exp[-2A;||^] (12) 
\ h 

where the sum is over all k» modes for which the AP test can be applied. From Fig. [I], 
we will take 0.002 < hi < 0.02 s/km, where the lower limit is set by the scale at which the 
slowly-fluctuating continuum would contaminate the signal, and the upper limit is set by 
the scale at which nonlinear distortions would start to become important (see [Hui 1998a] ). 



One can replace the summation by an integration: J2k n ~^ (L/tt) J dk\\ where L is the 
length of the absorption spectrum available. Assuming full coverage between Lya and Ly/3, 
L ~ 50000 km/s. 

Hence, taking model A to be the Q m = 0.3 - Qk = 0.7 universe and model B to be 
the Q m = 0.3 - = 0.7 universe, we find that, at z — 3, the expected S/N is 0.8 and 
0.6 respectively for angular separations of 1' and 2'. It turns out S/N peaks at about 0.5', 
reaching 1, but drops off at smaller #'s. (The S/N estimate above should be modified for 
sufficiently small #'s because P» 3> P^ no longer holds. Also, it changes somewhat with z, 
n and the fen's we include, but it provides a good rough estimate for a reasonable range of 
parameters.) There are two competing effects: a larger 9 gives a larger difference between 
models (Fig. |1]), but also a smaller P^, hence harder to measure. Note that this is for only 
one pair of quasar spectra. To reach a 4 — a level discrimination, something like 25 pairs, 
at angular separations 0.5' — 2', would be required. Note also that in the above estimate, 
we have ignored the error in the prediction of the cross-spectrum from the observed auto- 
spectrum. This error is likely to be sub-dominant, particularly because a very large number 
of LOS is in principle available for measuring the auto-spectrum. There are roughly 10 pairs 
of quasars with existing spectra at the above angular separations, or slightly larger, for z > 1 
(see e.g. |Crotts fc Fang 1998| & ref. therein). Upcoming surveys such as the AAT 2dF and 



SDSS are expected to increase this number by at least an order of magnitude. Some of these 
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will be faint quasars for which shot-noise might be important. 

A few issues should be further explored in the application of the AP test. First, it is 
obvious from the above analysis that we could boost the S/N for the smaller angular separa- 
tions by extending to higher fen's. This requires, however, the modeling of nonlinear redshift 
distortions. The finger-of-God on small scales would lower our predicted cross-spectrum. The 
exponential factor in eq. (|J) could be used for modeling this (e.g. Peacock fc Dodds 1994| ), 



but it probably depends on more than simply the £Ts we are interested in. Moreover, our 
large-scale linear distortion kernel (eq. ||) should be checked against simulations for accu- 
racy. It is known in the case of galaxy-surveys, for instance, that the linear prediction for 
redshift distortions could be modified even on very large scales. Furthermore, the effective 
bias factor in the distortion kernel (2 — 0.7(7 — 1) i n Pf] eq. |7]) might also deviate from the 
linear prediction ( |Hui 1555a ). 



Second, we have focused on a particular version of the AP test here, in which we 
construct a predicted cross-spectrum based on the observed auto-spectrum. In practice, one 
should test the whole inversion procedure from the auto-spectrum to the cross-spectrum 
with simulated noisy data, to guard against possible instabilities. An alternative would be, 
for instance, to parameterize pf(k) (eq. || & flTD| ) by an amplitude and a slope, and then 



fit simultaneously these parameters and the various O's to match the observed auto- and 
cross-spectra. This could lead to smaller error-bars for the measured JTs by restricting the 
form of Pf(k). 

Lastly, we have chosen to focus on a subset of cosmological parameters which are deemed 
to be currently popular (Fig. |2|). However, there exist cosmological models which are less 
conventional, but nonetheless not necessarily ruled out by current observations, such as a 
universe closed by a large Qq (e.g. |Kamionkowski fc Toumbas 1996| ). Such models could 
predict a cross-spectrum so strong that even one pair of quasar spectra would be sufficient 
to rule them out, that is, if they do not describe the actual universe. We will pursue this 
and other observational issues in a separate paper. 

As this work was nearing completion, we became aware of efforts by several groups who 
considered similar ideas flGroft 1998fc McDonald fc Miralda-Escude 19981 ; |Seljak 1998Q . This 



work was supported by the DOE and the NASA grant NAG 5-7092 at Fermilab. 
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Fig. 1. — The fractional error in the predicted cross-spectra (each panel for the cross- 
spectrum of the given angular separation) for 5 different assumed cosmologies, labeled ac- 
cording to the order of the curves from top to bottom. The input cosmology coincides with 
the assumed model of the solid line. The fractional error is defined as 5P£/ , where 

input 

5Pl = Pl-Pl. and P? is the predicted/output cross-spectrum. 7 is assumed to be the 

input 

same as the input value (see text). 



13 




Fig. 2. — The velocity-angle ratio (eq. 0) at a redshift of function of two 

parameters. The bottom solid curve shows u±/8 as a function of w (bottom axis) for models 
in which Vt m = 0.3 & Qq = 0.7. The upper dotted curve shows u±/9 as a function of fl m 
(top axis) for models with no Q, i.e. Q m + = 1. Models with mixture of Qq and flk 
smaller than 0.7 generally have u±/8 somewhere between the dotted and solid lines. 



